Phase diagram and critical exponents of a dissipative Ising spin chain 

in a transverse magnetic field 
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We consider a one-dimensional Ising model in a transverse magnetic field coupled to a dissipative 
heat bath. The phase diagram and the critical exponents are determined from extensive Monte 
Carlo simulations. It is shown that the character of the quantum phase transition is radically 
altered from the corresponding non-dissipative model and the double-well coupled to a dissipative 
heat bath with linear friction. Spatial couplings and the dissipative dynamics combine to form a 
new quantum criticality which is independent of dissipation strength. 



Dissipation in a quantum system is an important sta- 
tistical mechanical problem having ramifications in sys- 
tems as diverse as a single spin coupled to an environ- 
ment 0] to limitations of quantum computation 0- Of 
particular importance are the connections among deco- 
herence, noise, dissipation, and the amount of coarse 
graining necessary for classical predictability Q. Dis- 
sipation can arise when we focus on some distinguished 
variables of a system, which are coupled to the ignored 
environment variables that are integrated out in a closed 
universe. Since dynamics and thermodynamics are in- 
timately intertwined in quantum mechanics, dissipation 
plays a very important role. 

Although a single two-state system coupled to dissipa- 
tion has been extensively discussed in the literature [J], 
and its practical implications are abundant, the case of 
infinitely many spatially coupled two-state systems has 
not been discussed. We shall show that spatial couplings 
and the dissipative dynamics combine to determine a new 
class of quantum critical points, which is different from 
the dissipative phase transition of the single two-state 
system j5f and the corresponding non-dissipative model. 
The nature of the quantum criticality does not depend 
on the value of the dissipation strength. This finding is in 
striking contrast to predictions for an array of resistively 
shunted Josephson junctions 0,0- 

The Hamiltonian for the Ising spin chain in a trans- 
verse field and coupled to a heat bath reads 



N, 



A, 



1=1 i=l 



i,k 



where a x and a z are the Pauli matrices. To each site, i, 
we have coupled a set of independent bosons {a^fc} (de- 
struction operators), of frequency Wi,k, representing the 
environment. The coupling strength, Cfe, to the coordi- 
nate of the fcth oscillator is chosen such that the spectral 



function J(u>) = Air C%6(w — uji.k) = 2irctuj, for to less 
than a cutoff frequency uj c , but vanishes otherwise. This 
spectral density defines linear dissipation whose strength 
is determined by the parameter a [8j- 

The bath degrees of freedom can be integrated out in 
a manner identical to the treatment of the spin-boson 
problem The only difference is to recognize that 

the transfer matrix of the two-dimensional classical Ising 
model is the one-dimensional Ising model in a transverse 
field (non-dissipative, of course) [lfj. The partition func- 
tion defined on a (1 + l)-dimensional lattice indexed by 
i and r is 



Z = Z a 



(2) 



where Zq is the free boson partition function and 

s = —K E E w.TMt+i.T - r m,TtM,T+i 



i=l t=1 



i=l T = l 



i=l t<t' v ISt 1 



I)' 



(3) 



The classical Ising variables are defined by {/ii. T = ±1}, 
and the total number of lattice sites in the spatial and 
imaginary time directions are N x and N Tl respectively. 
Periodic boundary conditions are applied both in the 
space and imaginary time directions. The long-ranged 
term containing the dimensionless variable a describes 
the Ohmic dissipation Here K = t c J and e~ 2r = 

tanh(T c A), where r r is the lattice spacing in the imagi- 
nary time direction [13 ■ The universal scaling properties 
should not depend explicitly on r c , and we can set it to 
unity without loss of generality. We shall investigate the 
critical behavior of this partition function as a function 
of the parameters K, T and a. 

We use Monte Carlo simulations to study the action 
in Eq. 0, employing a variant of the Swendsen-Wang 
cluster algorithm [DJ developed by Luijten and Blote 
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FIG. 1: Data collapse of Binder cumulants at the critical point 
for a — 0.6 and z — 2. Points to the left of the dashed line 
were not considered since finite size effects become important 
in this region. 

|l4|. This cluster algorithm greatly reduces autocorre- 
lation times and allows us to simulate lattices of more 
than 10 7 spins. First, assuming a continuous phase tran- 
sition, we attempt the scaling procedure of Ref. (lif . 
For each value of a we determine the critical coupling K c 
and the dynamical critical exponent z self-consistently by 
looking for a data collapse of the Binder cumulant ratio 
B = 1 — (m 4 )/(3(TO 2 ) 2 ). Near the critical point B scales 



B(N x ,N T )=cj> 1 



Z (N x y 



(4) 



In particular, at the critical point, where the correlation 
length £ diverges, the Binder cumulants collapse onto a 
universal function of N T /(N X ) Z . The strategy is thus to 
plot B as a function of N T for different values of N x . At 
the critical coupling K c the maximum value of B(N X , N T ) 
is independent of N x , and the critical exponent z can be 
found from the optimal collapse of these data sets onto 
the universal curve <t>Fi(0, N T /N Z ), as shown in Fig. ^ 
We use re-weighting techniques |l8| to vary the value of 
K continuously and omit small system sizes where finite 
size corrections become important. The excellent data 
collapse validates the scaling Ansatz in Eq. Q. 

At first, we fix T = -\ ln[tanh(l)] = 0.136 and cal- 
culate the phase diagram as a function of a and K . It 
is shown in Fig. [3 and the corresponding values of K c 
and z are listed in Tab. [i] The critical point a* at 
K = was computed using the method in Ref. |23j to 
be a* = 1.2492(4). The phase boundary K c (a) seems to 
approach the K — axis with zero slope which points 
to an essential singularity. This can be understood from 
a simple mean-field argument, since the divergence of 
the susceptibility x(a — * a*_) at K = is exceedingly 
strong. The direction of approach to the critical point in 
the fugacity-a plane does determine the precise form [5|, 
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FIG. 2: Phase diagram of the dissipative quantum Ising chain 
in the space of nearest neighbor coupling K and dissipation 
strength ot for Y — 0.136. The phase boundary appears to 
approach the K = axis with zero slope, well approximated 
by an essential singularity. 



but the divergence has an essential singularity as a — > a*_ 
(X(a) - e A/V<x*-ay Therefore, for finite but small K c , 
the phase boundary of the coupled system is given by 
2x(a)A" c = 1, hence the result shown in Fig. El 

We next consider the order parameter correlation func- 
tion C(x,t), which in an infinite system scales as 



C(x,t) = (a XjT aofi) — (<t x ,t) {&o,o) 
- x- {z+ ^g{x/^T/x z ). 



(5) 



The spatial correlation length £ diverges at the critical 
point and the characteristic frequency scale vanishes as 
£~ 2 , as long as the dynamical critical exponent z is a 
finite and well defined quantity. At the critical coupling 
K c , Eq. (J5J simplifies to 

C(x,t) - x- {z+r >-^g(T/x z ) = T-t* +n - 1 V*g(x*/T). (6) 

If we set r = and K = K c , the spin-spin correlation 
function in the spatial direction should therefore decay 
asymptotically as C(x, 0) ~ x~( z+v ~ 1 \ Observing that 
N T > 5N Z is sufficient to avoid finite size effects in the r- 
directions we chose N x = 140 and N T = 10 5 and show the 
result in Figure Fits were performed to the functional 
form (x~ k + (N X — x)~ k ) in the region x m - m < r < N x /2 — 
70. The optimal value for k — z + r/ — 1 decreases from 
1.003 for x m in = 5 to 0.99 for x m i n — 20 while the error- 
bars increase in size. We estimate z + r\ = 2.00(1). 

The imaginary time correlation function, shown in 
Fig. 0] falls to zero exponentially on the ordered side and 
(cz.tCo.o) approaches (cr X T ) 2 . On the disordered side, 
the correlation function cannot decay faster than t~ 2 0] 
and does not allow one to define a correlation time and 
a dynamical critical exponent z except in the sense of 
a generalized "Josephson length" |2^|. Taking into ac- 
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TABLE I: Critical coupling and the critical exponents z and 
v for different values of the dissipation strength a and the 
coupling r. 



r 


Q 


K c 


z 


V 


0.136 





1 


1 


1 


0.136 


0.2 


0.457699(5) 


1.97(3) 


0.637(7) 


0.136 


0.4 


0.23396(1) 






0.136 


0.6 


0.103842(2) 


2.00(2) 


0.639(3) 


0.136 


0.8 


0.03142(1) 
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0.0031045(1) 


1.97(4) 
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FIG. 3: Decay of the spin-spin correlation {<J x , T a x + rtT ) at the 
critical point for a = 0.6. The size of the lattice is N x = 140 
and N T = 10 s . The curve shows the best fit to the data points 
in the interval 10 < r < 70. 
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FIG. 5: Plot of 21n[m 2 ] - ln[m 4 ] as a function of lnN x at 
a = 0.6 and fixing N T = 5(N X ) Z ). The quantity [m n ] is 
defined in Eq. JHJ and according to Eq. JJJ the slope of the 
fitted line gives \ jv. 



FIG. 4: Spin-spin correlation function in imaginary time for 
a = 0.6, N x = 500 and 7V T = 10 4 . The three curves show 
(from top to bottom) the decay at the critical point, and in 
the disordered and ordered phase respectively. 



count the observed large corrections to scaling, our re- 
sults at the critical point K = K c , shown in Fig. Q]for 
N x = 500 and N T — 10 4 are roughly consistent with Eq. 
(|fj|) and the previously determined values of z and r/. On 
the other hand, it was shown that for the corresponding 
(0 + l)-dimensional problem C(r) ~ 1/lnr at criticality 
[17|. Because of smaller corrections to scaling, we used 
the spatial dependence of the correlation function for the 
quantitative analysis. 

We finally determine the critical exponent v. By choos- 
ing N T = b{N x ) z we again reduce the two-parameter 
scaling to a one-parameter one and can employ the same 
procedure as for the classical Ising model 19) . From the 
scaling ansatz (m n ) = 3 /l '(f> n ({K - K c )Nl /v ) it fol- 



lows that 



2 In 



dK 



- In 



dK 



— const 



-hxN x . (7) 
v 



The derivatives of (m™) with respect to K can be calcu- 
lated from the correlation function 



d(m n ) 



dK 



= (m n E x ) Kc - {m n ) K (H x ) Kc =: [m n ], (8) 



where S x denotes the quantity Y. x = J2 X t ^^^i+i.t- 
Results for u, obtained from fits such as shown in Fig. 
arc listed in Tab. D The errors due to the uncertainty 
in z were estimated from simulations with system sizes 
N T = 5(N X ) Z , using both z = 1.98 and z = 2. Within 
error bars there is no dependence of the critical exponents 
on the dissipation strength. 

To test if our choice of the third parameter in the ac- 
tion in Eq. ® - the nearest neighbor coupling V - affects 
the values of the critical exponents, we repeated the sim- 
ulations for r = — | ln[tanh(0.1)] at a = 0.6 and found 
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TABLE II: Comparison of exponents obtained by our sim- 
ulations and the values predicted by the e-expansion. The 
values of v and z are averages of the various estimates listed 
in Tab. Q The sum z + r\ was determined for a = 0.6 only. 





simulation 


one loop two loop 




V 


0.638 ± 0.003 


0.583 [20] 0.633 ± 0(e 3 ) 


w 


n 


0.015 ±0.020 


0.020 ± 0(e 3 ) 


[20] 


z 


1.985 ±0.015 


2 1.980 + 0(e 3 ) 


[20] 


Z + T] 


2.00 ±0.01 


2 2 [20J 





consistency with the previous results. Note that when 
the Hamiltonian of the Ising model in a transverse field 
contains the coupling to the heat bath, T is an indepen- 
dent parameter. We finally checked that our numerical 
results do not depend on the aspect ratio by increasing 
N T /(N X ) Z from 5 to 20 and used two different random 
number generators [2^1 l2fi| to rule out the only remain- 
ing possible sources of systematic errors. 

We now compare these result with the predictions from 
a dissipative n-component cf> 4 field theory ji^, in an 
expansion in e = 2 — d (the one dimensional Ising chain 
corresponds to n = 1 and e = 1), which is based on their 
result that the dissipation strength does not renormalize 
to all orders in e. The agreement with the simulation 
shown in Table ILTI is very good, but may deteriorate with 
higher order terms, as even the e-expansion of the stan- 
dard (f> 4 theory is asymptotic in nature and reliable expo- 
nents cannot be obtained without a Borcl-Pade analysis 

M 

We believe that for the first time we have been able to 
carry out a precise simulation of an extended dissipative 
quantum system and have made progress towards under- 
standing of quantum criticality in the presence of dissipa- 
tion. An intriguing problem that can be analyzed by the 
present method is the criticality of a (0 ± l)-dimensional 
transition embedded in the higher dimensional system 
discussed in the context of an array of resistively shunted 
Josephson junctions 0,0 and the quantum theory of the 
smectic metal state in stripe phases |2^. It is remark- 
able, however, that a floating phase arises naturally in 
the problem of a dissipative Josephson junction array, 
where the states of (0 + l)-dimcnsional power law phases 
can slide past one another. 
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